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I. INTRODUCTION 


The purpose of this thesis is to examine two particular 
aspects of underwater acoustic scatter from the Arctic ice 
canopy. The motivation for this study came from previous 
laboratory experiments during which some unexplained 
phenomena were observed. 

The method used in this research is to supplement 
available theory of scatter by the use of scale models 
constructed of a material that effectively represents the 
known physical properties of Arctic ice. Understanding 
sound propagation in the Arctic is important for naval 
applications. In a recent Naval Postgraduate School thesis 
by Denny and Johnson (1986), a study of broad scope was made 
to determine the importance of various characteristics of 
the ice canopy on the overall level of low frequency 
backscatter from the ice. In the case of ice ridges and the 
ice edge at an open lead (or polynyas), the Biot-Tolstoy 
theory for diffraction from an infinite, rigid wedge was 
applied in an attempt to predict the amount of energy which 
would be backscattered from these acoustically significant 
features of the ice. 

The first series of experiments in the present work 
explore the sound backscattered from the edge of an ice floe 
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at a polynya. An open lead or polynya in the ice cover is 
an area where there iS no ice. Although definitions vary, 
leads are generally cracks in the ice of length much greater 
than the width of the crack, and polynyas are vast areas of 
open water beyond the edge of an ice floe. 

Experiments performed by Denny and Johnson (1986), and 
independent measurements performed for the purpose of this 
thesis, have shown that more energy than predicted by the 
Biot-Tolstoy theory may be backscattered from the ice edge 
of the laboratory model of a polynya. Evidence will be 
presented to show Arctic ice and the acrylic ice model are 
acoustically non-rigid, so some of the sound energy incident 
on the plate generates solid-borne compressional and shear 
waves which propagate along the plate and radiate energy 
back into the water when they reach a discontinuity, such as 
the end of an ice floe at a polynya. This reradiated 
energy, when superimposed with the actual diffraction, is 
the cause of the received backscattered signals being larger 
than that predicted by theory alone. 

The second set of experiments considers the sound 
backscattered from acrylic models of an ice ridge. The 
interaction of adjacent ice floes forms pressure ridges in 
the ice canopy. The thickness of the interacting ice and 
the magnitude of the compressive and shear forces exerted on 


the ice floes will determine the size of the keel depth 
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(below the ice canopy) and the sail height (above the ice 
canopy). For the laboratory model, the ice ridge is assumed 
to be a wedge of triangular cross section mounted below an 
acrylic plate. The sound energy backscattered from a wedge, 
as described by the Biot-Tolstoy theory, is the superposi- 
tion of the reradiated energy from each line along the crest 
of the wedge. The experiments were designed to determine 
what part of the wedge contributes the predominant backscat- 
tered energy, how the size and orientation of the wedge 
affects the amount of energy which is reradiated, and how 
reflection, diffraction, and mode-converted energy share in 
the total backscatter over a range of sound frequencies and 


incident angles. 
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Ii. RESEARCH FACILITIES AND EQUIPMENT 


A. OCEAN ACOUSTICS WAVE FACILITY TANK AND LABORATORY 

The Ocean Acoustics Wave Facility at the Naval 
Postgraduate School, previously used in the study of 
air-water interactions on the ocean surface, contains an 
anechoic tank in which this experiment was conducted. The 
tank was large enough to facilitate the use of sinusoidal 
pulses of energy to measure diffraction without interference 
from reflections due to the sides or bottom of the tank. 
The dimensions of the tank are 10 feet in length and width, 
9.5 feet in depth. The bottom of the tank and all four 
Sides are lined with redwood absorbing wedges which result 
in a low noise environment and less reverberation of sound 
in the tank. Signal averaging was performed using a very 
Short pulse repetition time in order to eliminate the 
effects of random background noise. 

The equipment used in the laboratory is listed in Table 
I. Figure 2.1 illustrates the equipment set-up used during 
the performance of both experiments. All of the equipment 
used is available on the commercial market with the excep- 
tion of the 12 element polyvinylidene fluoride (PVDF) planar 
transducer which was designed for future waveguide experi- 
ments, but found some application here. The design 
and operation for the PVDF transducer is described in 
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Abbreviation 


Timing Simulator 


Function Generator 


Power Amplifier 


Source or Receiver 


Pre-Amplifier 


Nicolet 


Voltage Amplifier 


Signal Processor 


Scope 


TABLE I 


BOUIPMENT Gist 


Nomenclature 

Interface Technology Timing 
Simulator/Word Generator, Model RS- 
648 


Hewlett Packard Function Generator, 
Model HP-3314A 


Hewlett Packard Power Amplifier, 
Model HP-467A 


Celesco LC-10 Hydrophone 


Ithaco Low Noise Pre-Amplifier, Model 
e20n: 


Nicolet Dual Channel FFT Analyzer, 
Model HP-465A 


Hewlett Packard Voltage Amplifier, 
Model HP-465A 


IBM PC/XT with Computerscope 
ISC-16 Data Acquisition and Analysis 
Package 


Kikusui Oscilloscope, Model COS 5060 
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HP-467A 
Pwr Amp 


$90006n 60 6us Gn COE EP OT 8 5599 0° aa eee ee eae mn) et eeN ee tne eS o” Soccer esseseoesccons 


ITHACO 1201 
Low Noise Pre-Amp 
OOo 
HP-3314A ° 


rong) Function Generator 
o 









==) 
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Se 
HP-455A 
Volt. amp 
00000000000C¢000 
SeEGoyooo 
Interface Technology 
Timing Simulator /Word Generator 
eaaeC 7 AT 
| 
Seo Oo 
G10 701-O 0 
Sono noTo 
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OG 0 10.70 





A/D Converter 











sampling frequency 
= 1 MHz 


Figure 2.1: Schematic diagram of the signal processing flow 
path. 
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Appendix A. Additional information concerning the equipment 
used for transmitting and receiving acoustic signals, and 
for performance of the signal processing are described in 


Sections B and C of this chapter. 


B. TRANSMIT AND RECEIVE SYSorteMs 
1. Hewlett Packard Function Generator, Model 3314A 

The 3314A function generator is a digital, multimode, 
programmable signal generator with preset Sine, Triangle, 
and Squarewave functions and programmable arbitrary 
waveforms. The frequency range is from 1 mHz to 19.99 MHz, 
and the maximum output voltage is 10 volts zero-to-peak. 
While using repeated bursts of sound energy for the purpose 
of signal averaging, the function generator was operated in 
the N CYCLE MODE, in which N complete cycles ofthe 
designated waveform were generated by the function generator 
at the start of each pulse repetition period. 

2. Nicolet FFT Wave Analyzer, Model 660 B 

The Nicolet wave analyzer is a self-contained, 
programmable signal processing system which performs analog 
to digital conversion internally. The display screen allows 
real-time observation of either the received time signal or 
its frequency spectrum. The Nicolet wave analyzer is 
limited to a sampling frequency of 250 kHz and a time window 
of 4 msec. When utilized in parallel with the IBM signal 
processing equipment, the Nicolet was useful for initial 


detection of the signal of interest. 
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Se Interface Technology Timing Simulator, Model RS-648 
The timing simulator was used to establish the pulse 
repetition rate used for signal averaging by triggering the 
function generator, the Nicolet, the Ithaco pre-amplifier, 
and the signal processing equipment in the proper sequence. 
The timing simulator has the ability to produce TTL pulses 
sufficient for a 50 nanosecond timing resolution. 
4. Ithaco Low Noise Pre-Amplifier, Model 1201 
The Ithaco pre-amplifier provided the filtering and 
amplification necessary to detect the small magnitude, 
coherent transient signals backscattered from different 
features of the acrylic ice model. Signal filtering in the 
bandpass range of 3 kHz to 100 kHz helped to reduce the 
level of background noise, and signal amplification at gains 
ranging from 100 to 10,000 times magnification permitted 
detection of even the smallest signals which were reradiated 
from the model. 
eee tesco EC—10 Hydropnone 
Celesco LC-10 hydrophones were used both as 
transmitters and receivers. The LC-10 hydrophone is small 
enough to simulate a point source of sound, in which ka << l 
(where k = wave number, a = linear dimension of the 
transducer). The LC-10 hydrophone has a fast response; 
transducer rise and decay times are short enough to permit 
detection of short pulsed signals free of transducer 


reverberation or ringing. Short duration pulses were used 
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to ensure that received signals from propagation paths of 
similar length would be separable in time and free of 
unwanted interference. The LC-10 also has the advantage of 
a flat free-field voltage response and an omnidirectional 
directivity pattern within +2 dB over a frequency range of 
30 kHz to 100 kHz. When measuring the relative amplitude of 
different signals at a given frequency, it was not 
necessary to know the absolute calibration of the 
hydrophone. The LC-10 does have low output power, but this 
disadvantage was not a limitation because the hydrophone was 


used at short ranges in a low ambient noise environment. 


C. SIGNAL PROCESSING EQUIPMENT 

The primary signal processing for the experiment was 
performed on an IBM PC/XT personal computer using a 
Computerscope Model ISC-16 Signal Processing system. The 
Computerscope ISC-16 system is a fully integrated hardware 
and software package designed to permit the IBM PC to 
perform as a data acquisition and analysis laboratory 
instrument. The ISC-16 system contains a 16 channel analog 
to digital converter, an external instrument interface, and 
the Scope Driver software. The system can receive 16 
channels of data input at an aggregate sampling rate of l 
MHz. Digital conversion with 12 bit accuracy is achieved 
over an input range of +10 volts. The Scope Driver software 


permits the IBM PC to operate like a digital storage 


ZZ 


oscilloscope with the ability to obtain graphical output on 


amedot Matrix printer or an X/Y plotter. 
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III. DESCRIPTION OF PPHYSICAL Merrie 


In order to accurately model the way sound in the Arctic 
Ocean interacts with the ice cover, it was necessary to know 
the physical properties of the ice and to find a material 
with similar properties to use in the laboratory. Since 
information available from the literature of Arctic 
Scientists is extensive (but sometimes conflicting since 
conditions in the Arctic are so variable) it was possible to 
compile data for the average material properties in the 
ArcGtre:? For a model designed to simulate the scattering 
properties of low frequency sound from the ice canopy, the 
typical values of the characteristic properties of sea ice 
are sufficient to acoustically model the backscattering from 
Arctic ice, since local variations in the ice properties 
typically extend over lengths small in comparison to a 
wavelength. 

Acrylic was chosen as the material for the model because 
most of its bulk properties fall within the range of 
physical properties measured in Arctic sea ice. The 
longitudinal wave speed and flexural wave speed (scaled for 
the appropriate thickness), the Poisson ratio, and the 
characteristic impedance (pc) for acrylic are all very 
Similar to ice. The acrylic plate density is about 25% too 
high, and compressional and shear wave attenuations in the 
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acrylic are higher than in ice, but these differences are 
acceptable. Table II summarizes the comparison between 
Arctic ice conditions and the laboratory model of an acrylic 


sheet floating in a fresh water tank. 


Pee MODELS OF SMOOTH ICE AND ARCTIC POLYNYAS. 

For a scale model to accurately represent the real 
physical situation, there must exist physical and geometric 
Similarity between the two situations. Physical similarity 
was assured when acrylic was chosen as the modeling 
material, since its bulk properties closely resemble those 
of Arctic ice, as described in the last section. Geometric 
Similarity is achieved by maintaining a constant ratio of 
wavelengths to thicknesses and distances between the Arctic 
and the laboratory situations. Medwin et al. (1984) 
successfully demonstrated the accuracy of this modeling 
technique, and similarly in this experiment, the ratio of 
the acoustic wavelength (i) to the ice thickness (h) must be 
the same for the Arctic and the laboratory situations. 

The low frequency range of interest for sound propaga- 
fon ih the Arctic ocean is 50 to 400 Hz. Assuming a 
typical Arctic level ice thickness of 2 m, then the 3/8" 
acrylic plate (9.5 mm) used as a model of smooth ice 


corresponds to a scale ratio of approximately 200:1. 
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Consequently, for the 9.5 mm plate, the following frequency 


matios apply: 


Arctic frequency (Hz) Lab frequency (kHz) 

pO HZ FO" kHz 

BOG Hz 80 kHz 
Heme the 178" (3.175 mm) acrylic plate, a similar corre- 
spondence results. Also, the scaled wavelength to thick- 


ness ratio will be met for the two situations since the 
speed of sound, c, in water agrees within 2% between Arctic 


ecean conditions and the fresh water in the laboratory. 


Then, 

wavelength (A) = speed of sound (c) / frequency (f) 
ene (arctic) ~ “(laboratory ) 
. BO aee 7 Beacon: 


Since the frequency and thickness were scaled by the same 
hactor. 

In the low frequency regime, if the roughness elements 
on the smooth Arctic ice are small compared to a wavelength, 
they are not acoustically significant and can be ignored. 
For the same reason, it was not necessary to introduce 
roughness elements on the acrylic model of smooth ice. The 
Arctic polynya, the large expanse of open water beyond the 
edge of an ice flow, was modeled by a single edge of a large 
acrylic plate. The length of the plate edge was 1.5 m (10 


to 80 acoustic wavelengths, depending on the frequency used) 


oil 


so that the ends of the plate would not interfere with the 
diffraction from the edge of interest. 

Finally, since the density of the acrylic is greater 
than the density of water, it was necessary to build-up the 
edge of the plate to keep it afloat. By gluing a 3/8" high 
lip around the outside edge of the acrylic sheet, the plate 
resembled a flat-bottomed boat floating on the water. Since 
the lip was glued to the top side of the plate, the edge of 
the plate which caused the incident sound to be diffracted 


was undisturbed. 


B. MODEL OF A PRESSURE RIDGE 

When pressure ridges are formed by the interaction of 
adjacent ice floes, the keel depth will naturally be a 
function of the ice thickness which formed the ridge. The 
data compiled by Hibler et al. (1972) from visual measure- 
ments above the ice, and the submarine underice profiles 
from below, resulted in a statistical model of pressure 
ridge keel depths and ice keel spatial density. To ensure 
that only pressure ridges are considered in the profiles, 
the data distribution only included pressure ridges which 
were greater than a predetermined cutoff height or depth. 
The average pressure ridge keel depth reported, from the set 
of all ridges greater than 6.1m, was 9.6 m in the Central 
Arctic and 11.4 m in the Canadian Archipelago. The average 
reported spatial distribution of pressure ridges was 4.3 
ridges/km in the Central Arctic and 9.6 ridges/km in the 


Sis 


Canadian Archipelago. Unfortunately, no data were reported 
for the average length of an ice keel. This type of data is 
not available from submarine underice profiles. 

For the purposes of this experiment, three different 
pressure ridges were built to study the effects of height 
and slope on the amount of energy backscattered from the 
ridge. Figure 3.1 shows a cross sectional view of the model 
for the "typical" pressure ridge, including the significant 
dimensions and the scaled dimension. It is the same ridge 
used by Denny and Johnson in their study of the backscatter 
from a ridge. The dimensional data for all these pressure 
ridges are presented in Table III. Wittheethe "ty pascal” 
pressure ridge as a reference, the "large" ridge was made 
with the same slope but twice the height, while the "steep" 
ridge was made with the same height but twice the slope. 

The scaled dimensions of the "typical" pressure ridge 
were designed to match those of a typical multi-year 
pressure ridge, while the "steep" and "large" ridges 
represent possible ridges specifically designed to study the 
dependence of backscatter on height and slope. Each keel 
was manufactured by gluing flat strips of acrylic together, 
then machining the entire piece to the specified dimensions. 
The ridges were glued to a smooth acrylic sheet using an 
acrylic solvent (K-Lux™ solvent cement for acrylic). 

It is noteworthy that the fabrication of the "steep" and 


the "large" ridges did not include a ridge sail above the 
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DIMENSION 


Keel Base 


Height 


Angle 


Sail Base 


Height 


Angle 


TABLE SET 1 


COMPARISON OF THE DIMENSIONS 
OF THE THREE PRESSURE RIDGE MObEGe 
FOR A 1/8" PLATE THICKNESS 


RIDGES 
TYPICAL LARGE STEEP 
RIDGE RLVGE RIDGE 


2-3/8" (36.2 m) 4-3/4" (72.4 m) 27/32" (18 23eeae 


5/ 8s (955) mm) S117 43 ieee 5/8%. "(9 soa) 
28° Zo 565 

3/45) (a4 me (N/A) (N/A) 

3716 ie (2 iy (N/A) (N/A) 
26 (N/A) (N/A) 


(Dimensions in parentheses are scaled to the Arctic model 


assuming 


a frequency ratio of 600:1.) 





Figure 3.1: Cross section of the pressure ridge model used 


to measure forward and back scatter from a 
ridge keel. Dimensions in parentheses are 
scaled to the Arctic model. All materials are 
aecry lier 
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plate. Although the sail is a visually predominant feature 
in the Arctic (and on the "typical" ridge), it is assumed to 
be acoustically insignificant for the purpose of diffraction 
and backscatter from the ice keel. Additionally, it was not 
necessary to place roughness elements on the faces of the 
pressure ridges since, for low frequency propagation of 
sound in the Arctic, the roughness elements are small 
compared to a wavelength and, therefore, are also assumed 
acoustically insignificant. This may not be the case and 


will be examined in future experiments. 
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IVse 0 HEORY 


A. DIFFRACTION  THGORyY 

A physical explanation for the mechanism which causes 
diffraction is best described using Huygens! principle, as 
described by Clay and Medwin (1977). Huygens' principle 
states that each point of an advancing wavefront is a source 
of secondary waves which move forward as spherical wavelets 
in an isotropic medium, and the outer surface enveloped by 
all of the wavelets defines the new wavefront. Consequently, 
when an acoustic wave interacts with an object in the 
medium, every point on the object will be the site of an 
expanding spherical Huygens' wavelet. The wave front that 
envelopes all of the wavelets spreading from the object is 
called the diffracted wave. 

Biot and Tolstoy (1957) derived a closed-form solution 
to predict the maximum value of the diffracted energy from 
an infinite rigid wedge due to an impulse point source of 
known strength. The expression was simplified by Medwin 
(1981) to consider a point source delta function turned on 
at time t = 0, and to account for the finite size of the 
wedge using digital techniques. The result is the 
Simplified expression in cylindrical co-ordinates for the 


diffracted pressure from a wedge, Equation 4.1: 
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p(t) = [-(Soc/478,,){8texp(-ty/4,,)]/(rro sinh y) fairl} 


(4,2) 
Suerte wine C. ) | 





{a} = 
ee ce xD eaamny / C...) cos(7/6,.)(™6+0+6,)]) + ExP(sZty/ Gs) 


2 2 2 
to = (r° + ee + 27) 


a sil 
y cosh ; ; (4.3) 
rro 


The diffracted pressure at time t is p(t), op is the density 





of the medium, c is the sound speed in the medium, S is the 
source strength. The angle of the wedge measured in the 
Pluic region is €,, the source coordinates are (ro, fo, 0) 
Eemomecne receiver coordinates are (r, %, Zz}. TRewmeer mn 2. Or 
Mem@erce written for Simplicity as (7+6+€,), and actually 
represents the sum cf the four possible combinations of +8 
and +*,. Figure 4.1 shows the geometry for the wedqe. 

For a point source turned on at time t = 0, the sound 
Giverqes spherically and intersects the ridge first at the 
least time path, then along the crest of the ridge as the 
wave diverges. Sound energy diffracts all along the crest 
of the ridge and propagates to the receiver according to 


meygens' principle. The earliest arrival occurs at the 


least time, To, given by 


Bs 


mae ae ZC 


Se 


(a) 





wt 


+ Ow 


a 
wevrvgy 
\ 


RCYR 
Least Time 
Z=0 Path 
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oe bd +o eee Ae 


fl 


0 ee oe 


cS 


Diffrecting 
Edge 


Piguce: seule a) Wedge aeometry for equations 4.1 through 


(a) 
4.3. (b) Unfolded geometry. 
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The maximum pressure is received at this least time, while 
pressure received from other paths falls off exponentially 
with time corresponding to the greater path length and 
weaker diffraction at oblique angles. 

Novarini and Medwin (1985) demonstrated the use of the 
"wedge assemblage" method as an accurate way of modeling the 
diffraction from a conjoint assemblage of rigid or pressure- 
release crests and troughs. Similarly, cine total response 
for the diffracted pressure from the three Biot-Tolstoy 
wedges which form the ice keel in the acrylic model will be 
the linear superposition of three properly timed wedge 
impulse responses (assuming for simplicity that the wedges 
are rigid or pressure-release). Figure 4.2 shows the 
pressure due to three rigid wedges which form the model of 
ine ice keel. The wedges are also assumed to be 
sufficiently spaced to permit superposition of the 
individual responses. 

In the case of backscatter, the geometry is simplified 
pce Fr — Yor 8 = 65, and z = 0. imetiesseecase , £ rOm 


Equation (4.3), 


Car) = 2r" 
y(t = tT 5) = cosh~) | ———————— | = cosh“'(1) = 0 
Ores 
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(a) 


WLLL LLL dd hdd, Ui SLOT TTT, 


x (3) 








YY 


(b) 


p(t) 





Poo 6 60° 120 180 240 Seo 
Wane (see) 


Fiqure 4.2: (a) Diagram showing the three Baot-Tolseua 
weddqes that make un the model of the ridge 
keel. (b) Time domain impulse response 


calculated using Biot-Tolstoy theory age 
superposition. Relative times of arrival fora 
backscattering qeometry. 
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pence spl(t) <{[sinh (y)J~*, plt) > = at t = to, aS expected. 
Since digital signal processing is used, Medwin (1982) has 
Shown that, in many cases, the initial impulse can be 
approximated by a large finite value corresponding to the 


adjacent sampling interval, At: 


p(t = To) 
= 1.366 
pitts tG teat) 


Banally, in spite of the fact that the Biot-Tolstoy 
theory predicts the diffracted energy from a wedge due to an 
impulse of sound, the theory can be used to predict the 
diffracted energy when the input is not an impulse by 
relying on linear systems theory. If the ice keel is 
considered to be a linear system whose output is the 
diffracted pressure, then for an impulse source 4(t), the 
diffracted pressure given by the Biot-Tolstoy theory is the 
impulse response h(t). The transfer function of the filter, 


mee, 1S the fourier transform of h(t). 
ot) ete 
ae }-—3 
x ce yGe) 
In the experiments, a pulsed sinusoid x(t) is used as the 
input instead of an impulse source, so the output diffracted 
pressure response is the convolution of x(t) with h(t). 


prey — x(t) * nit) 


41 


Converting to the frequency @domatn oy the convolution 
theorem: 


P(f) = (X(£) aE Ce 


As shown by Denny and Johnson (1986), the pressure 
response in the frequency domain relative to the direct 
free-field transmission at a reference range R is the filter 


transfer function times the reference range: 


Po(f) XC fe} 
Pe Si) er 


REL 
Peet) X(£)/Ro 





= Fiore 


Consequently, for a known input Signal to the linear filter 
H(f), the output response can be obtained relative to the 
response over a reference distance even if the input was not 


an impulse. 


B. MODE CONVERS ONS aL ory 

When a sound wave interacts with a solid structure, some 
of the energy may be transmitted from one medium to the 
other causing longitudinal and transverse vibrations in the 
solid, and this process is sometimes caiied mode conversion. 
As described by Beranek (1971), when steady state is 
reached, the vibration field will build up to a level at 
which the power from the source balances the power losses in 
the system. If a burst of sound interacts with a solid, a 


Similar energy balance can be performed to account for the 
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disposition of the incident energy which gets reflected, 
transmitted, or causes vibration in the solid. 

For a wave to propagate in a solid, liquid, or a gas, 
the medium must be capable of alternately storing energy in 
kinetic and potential forms. Potential energy is stored 
when part of the medium has undergone elastic deformation; 
kinetic energy is stored in the movement of particles which 
have mass. A fluid does not support shear motion, but can 
undergo compression, so only compressional waves exist ina 
Paid . 

A solid can store energy in both compression and shear, 
so many different types of waves (representing various 
combinations of compression and shear) can propagate in a 
solid. For the flat plate used in the model of the Arctic 
ice, the solid or body waves of interest are compressional, 
shear, and flexural waves. A compressional wave is a 
longitudinal wave in the solid, such that deformations occur 
parallel to the direction of propagation of the wave, and 
the potential energy is only stored in compression. A shear 
wave 1S a transverse wave in the solid, such that 
deformations occur perpendicular to the direction of 
propagation of the wave, and the potential energy is only 
stored in shear. A flexural wave is a bending wave in which 
the potential energy is stored in a combination of 
compression and shear as the plate bends about its neutral 


axis (where there is no deformation). 
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1. Compressional Waves 


When a compressional wave propagates ina solid body 
whose dimensions are large compared to a wavelength, the 
wave will travel at the bulk speed of sound which is 


dependent only on the material properties of that solid: 





hy 
Y(l - o) 
Ca = 
6 (10) (eae) 
where Y = Young's Modulus (N/m*°), 6 = Poisson's Yatlonmomee 


density of the material (kq/m’). For the acrylic, vusamemene 
values from Table II, cp = 2545 m/sec. By Snell's Law, an 
incident plane wave at the boundary of the ice-water 
interface will be refracted into the second medium according 
to the ratio cf the sound speeds. Assuming that Huycens' 
principle applies, the wavelet initially produced at the 
acrylic boundary from the incident plane wave in the water 
responds to the bulk speed in the acrylic since it is 
unaware of the finite plate thickness, as shown in Figure 


4.3. In this case, the critical ang toe, sar 


Ss Cw _ 37° 
ae = sin?! —— = s1n a ee 
Cc, CB 


(grazina angle = 53°) 


At this critical angle, the energy transmitted into the 
acrylic is parallel to the surface, and will propagate as a 


longitudinal wave at the plate speed given by: 
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Cp — a ee oe en) Sec 


Similar results occur in an Arctic ice sheet at angles and 
speeds which maybe determined for the wide range of physical 


parameters that characterize Arctic ice (see Table II). 


8 
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Figure 4.3: Transmission from one medium to another at 
oblique incidence. 


Coe CO ne ec Ounce hn water 

Cp = compressional wave speed in a plate 
Cs = shear wave speed in a plate 

8; = angle of incidence 

8+ = angle of transmission 


(NOTE: Reflected rays, which occur in all cases, are not 
shown for simplicity.) 
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2. Shear Waves 

A shear wave will be generated in a solid in the 
same manner described above for compressional waves if the 
shear speed is greater than the speed of sound in the 
adjacent medium, which is the case for shear waves in much 
of the Arctic ice. In the acrylic model, the shear wave 
speed is less than the speed of sound in water, so 
refraction into the plate occurs at angles less than the 
incident angle. (Figure 4.3) Since a shear wave iS a 
transverse wave, optimum excitation of shear waves will 
occur when the acoustic energy in the water interacts with 
the floating ice plate or acrylic sheet at normal incidence. 

Considering Huygens' principle, for a plane wave in 
water at normal incidence, the wavelets transmitted into the 
acrylic cause vibrations perpendicular to the plate surface, 
resulting in transverse waves that propagate along the plate 
at the shear wave speed, cg = 1373 m/sec. In the Arectiter 
the typical shear wave speed is greater than or equal to the 
speed of sound in water, so shear waves in the ice plate 
will be excited by the same method described for 
compressional waves. 

3. Flexural Waves 

The theory of flexural waves and the derivation of 

flexural wave phase and group speeds are presented in 


Appendix B. A few points are in order here pertinent to the 


study of mode conversion. 
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When an incident acoustic wave strikes the plate 
causing vibrations, flexural waves may begin to propagate, 
in addition to the waves described above. The boundary 
condition for sound transmission between two media is that 
the horizontal wave numbers, k¢ in the plate, and ko in the 
fluid, must be equal to satisfy Snell's Law. When k¢ = ko, 
the compressional wave sound speed in water equals the 
flexural wave speed in the ice, and the situation is called 
coincidence. In the case of an acrylic plate or Arctic ice 
floating on water, the flexural wave speed is always less 
than the speed of sound in water, so the coincidence 
condition is never satisfied. Consequently, flexural waves 
cannot be efficiently generated by an incident acoustic 
plane wave. However, a source of diverging waves near a 
plate containing surface roughness elements can generate 
flexural waves (Denny and Johnson, 1986). As an additional 
consequence, if a flexural wave is generated in an acrylic 
or Arctic ice plate, it cannot reradiate energy efficiently 
back into the water in the form of compressional waves 
(except perhaps at a discontinuity) since the horizontal 


wave numbers cannot be matched. 
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V. FLEXURAL WAVE ATTENUATION 


The ineffective transfer of acoustic energy from the 
waterborne compressional wave to the ice medium in the form 
of plate flexural waves was studied by Denny and Johnson 
(1986). They showed that flexural waves could be driven in 
the plate and could radiate energy back into the water only 
when the source was very close to the plate. The flexural 
wave speed was effectively measured in that experiment, but 
no attenuation measurements were possible. The purpose of 
this experiment was to measure the flexural wave 


attenuation. 


A. EXPERIMENTAL PROCEDURE 

These measurements were made with the smooth floating 
plate in the anechoic tank to simulate the level ice model. 
In the flexural wave speed measurements made by Denny and 
Johnson (1986), an accelerometer mounted on a small square 
of plexiglass and bonded to the plate using Tackiwax™ 
provided an effective means of measuring flexural wave 
speed. While using the accelerometer, the plate-to-receiver 
efficiency depended on the quality of the bond. In order to 
measure attenuation, the receiver had to be constantly 
repositioned to find the decay of the propagating flexural 
wave over increasing distances. Since the method of uSing 
Tackiwax to stick the accelerometer to the plate could not 
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be duplicated from range to range with a measurable 
precision, the accelerometer could no longer be used as a 
receiver. 

The receiver used for the attenuation measurements was a 
Bruel and Kjaer 1/2" condenser microphone. The source used 
to excite the flexural waves was a 1" cylindrical hydrophone 
taken from a sonobuoy, and chosen for use because of its 
flat face which provided an efficient bond with the 
underside of the plate. The electrical input to the source 
was 3 cycles of a 200 volt peak-to-peak sinusoid. The 
frequency range was limited between 10 kHz and 20 kHz: the 
lower limit was determined by the maximum allowable pulse 
length for separating signals in time, while the upper 
frequency was limited by the microphone frequency response 
moelk-offt. 

The microphone was suspended from a metal rod spanning 
the plate and was slid across the plate at a height of 2 mm. 
Three data runs were conducted, one each at 10, 15, and 20 
kHZ. For the data collected at 10 and 15 Khz, the receiver 
waS positioned at 10 cm intervals between 1.0 and 1.5 m from 
the source. At receiver ranges less than 1.0 m, the time of 
arrival of the airborne sound wave generated by the source 
(travelling at 343 m/sec) interfered with the plate flexural 
wave (travelling at speeds less than 650 m/sec) for 
frequencies less than 15 kHz. At 20 kHz, the attenuation 


was too great for measurements to begin at 1.0 m from the 
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source, but the faster flexural wave eri (750 m/sec) 
provided sufficient signal separation to collect data 
between 0.5 to 1.0 m. Since the measurements were made in 
air, it was necessary to collect the data at night when the 
ambient noise conditions were low. The data were collected 
and processed using a Hewlett Packard 3561A Dynamic Signal 
Analyzer, by rms time averaging 1000 sample points to 


determine the received signal level. 


B. RESULTS AND Y©iSCUSsien 

As revealed by the experiments of Denny and Johnson 
(1986), the condition for efficient mode conversion between 
compressional waves in the water and flexural waves in the 
acrylic plate (or Arctic ice) is never satisfied. Bom 
efficient mode conversion, the horizontal component of the 
wave number must be matched between the plate and the water, 
but the condition cannot be met since the flexural wave 
speed is significantly lower than the speed of sound in 
water at all frequencies. Nevertheless, by using the 
transverse accelerometer, it is possible to excite flexural 
waves in the plate, albeit not efficiently. 

Figure 5.1 shows the measured flexural wave group speed 
based on the pulse time of arrival for the three frequencies 
tested. The measured wave speeds are less than the 
theoretical flexural wave group speed due to the mass 
loading effects of the water on the plate, which is 
unaccounted for on the theoretical curve. 
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magure 5.1: Flexural wave speed measured in a 1/8" 


acrylic plate floating on water. 
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The attenuation characteristic of the plexiglass was 
determined by measuring the largest peak-to-peak amplitude 
of the flexural wave signal, which corresponded to the 
maximum strength of the flexural wave at that frequency and 
range from the source. The attenuation (dB/m) for a 
particular frequency was determined from semi-log plots of 
received magnitudes versus range, such as in Figure 5.2. 
Data at longer ranges (e.g., 1.3 m in Figure 5.2) were 
affected by the noise level and were not considered in the 
determination of attenuation values. The data obtained at 
10, 15, and 20 kHz are plotted on Figures 5.2, 5.3, andmoe 
The cumulative results of those three figures form the 
attenuation versus frequency curve, Figure 5.5. 

The average attenuation of the flexural wave in the 
acrylic plate was measured to be 1.5 dB/m times the 
frequency in kHz. Comparison with sea ice values could not 
be made as no attenuation values could be found in the 


available literature for Arctic ice flexural waves. 
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Figure 5.5: Attenuation as a function of frequency for 
flexural waves ina 1/8" acrylic sheet floating 
on water. 
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VI. BACKSCATTER AT A RIGHT ANGLE WEDGE 
(ARCTIC OPEN LEAD) 
A. BACKSCATTER FROM AN UNDERWATER WEDGE 
The objective of this experiment was to measure the 
backscatter from the corner of an acoustically non-rigid, 
submerged, finite plexiglass wedge, and compare the results 
weem the predictions of the Biot-Tolstoy impulse solution 
moeemcdtitraction of point radiation by an infinite, rigid 
wedge. Measurements of backscatter from a floating acrylic 
wedge made by Denny and Johnson (1986) indicated that the 
Biot-Tolstoy diffraction theory does not accurately predict 
the amount of backscattered energy. In fact, in some 
directions more energy was received than predicted, but due 
EO the fact that the acrylic was floating, the effects of 
the air/water interface near the diffracting corner were 
unclear. By completely submerging a plexiglass wedge and 
repeating the measurements, the interactions with the water 
surface were eliminated so diffraction alone could be 
measured. 
mea xperimental Procedure 
The laboratory measurements of the backscatter from 
the corner of the wedge were performed in the Ocean Acous- 


tics Wave Facility tank described in Chapter II. The wedge 


ay 


was fabricated out of two pieces of 1/8" Sleminee (70 x 362 
cm). Liquid acrylic plastic cement was used to bond the 
corner of the wedge at a 90° angle; support pieces were 
cemented to the top and bottom (Figure 6.1). The wedge was 
suspended vertically in the anechoic tank by four ropes 
above and was held stationary by two side strings to prevent 
EOtat Tone, 

The LC-10 source and receiver were attached to an 
aluminum rod which was held vertically near the wedge by an 
extending arm from the side of the tank. Positioning of the 
source and receiver was accomplished using a 25 cm radius 
compass marked in 15 increments and centered above the 
corner of the wedge. The geometry for source/receiver 
positions was chosen to eliminate interference from 
undesired edges, since the objective was to measure diffrac- 
tion only from the wedge corner. The input electrical 
Signal used was 2 cycles of either a 90 or 100 kHz sinusoid. 
The received diffraction was identified both by its time of 
arrival and by disturbing the diffracted signal by means of 
an aluminum sphere placed at the edge. Figure 6.2 shows a 
typical received waveform, including the direct source-to- 
receiver path, the specular reflection from the plexiglass 
plate, and the backscattered signal from the corner of the 
wedge. Figure 6.3 represents the same geometry with the 
exception that an aluminum ball has been placed at the edge, 
and the diffracted signal can be identified by the 
additional response preceding the diffracted wedge signal. 
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Figure 6.1: Experimental set-up for the underwater wedge. 
Source /Receiver pestti190ned for monostatic 
backscatter, (r = Fo), (8 = 86,4). Angle of 
wedge measured in the fluid region, 4, = 270°. 


5 


Figure 6.2: 
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Backscattered waveform from an underwater 
wedge. 
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Backscattered waveform from an underwater 
wedge. 


1) = Direct path between source and receiver 

(2) = Specular reflection from the side of the 
wedge 

(3) = Diffraction from the edge of interest 

(4) = Diffraction distorted by the aluminum 


sphere 
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2. Results and Discussion 

Once the backscattered signal was identified, its 
measurement was simply the maximum peak-to-peak amplitude of 
the received signal. Backscatter wasS measured at angles 
me og = 15°, 30°, 45°, 105°, 120°, 135°. It was impossible 
to measure the diffraction between 50° and 95° because the 
desired signal was hidden by the stronger specular reflec- 
tion from the plexiglass plate. It was unnecessary to 
measure the backscatter for angles greater than 135° because 
of the symmetric geometry. The backscattered signal was 
identified using the aluminum sphere for each repositioning 
of the source/receiver. 

The results of the measured backscatter are plotted 
in terms of diffraction loss in dB referenced to spherical 
divergence over the same (r + rg) range as a function of 
angle from the side of the wedge (Figure 6.4). A single 
theoretical curve for the Biot-Tolstoy predicted diffraction 
loss from a rigid wedge is plotted for comparison. (A 
Single theoretical curve is sufficient because the predicted 
diffraction loss at 90 and 100 kHz are almost identical). 
Notice that the measured backscattered amplitude is greater 
than predicted for angles corresponding to backscatter from 
the trailing edge of an ice lead (60 = 6, = 15°, 30°, 45°). 
Also note that the measured signal is less than predicted 
for those angles corresponding to large backscatter from the 
leading edge of an ice lead (8 = 8, = 105°, 120°, 135°). 
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The experimental results shown in Figure 6.4 can be 
explained by the penetrability of the wedge causing mode con- 
version. For the angles corresponding to backscatter from the 
trailing edge of an ice lead (8 = 6, = 15°, 30°, 45°), some of 
the acoustic energy incident on the plate face is mode- 
converted and causes both longitudinal and transverse 
vibrations in the plexiglass. Compressional and shear waves 
propagate in the plate and radiate energy back into the water 
upon reaching the corner of the wedge. The speeds of propaga- 
tion are such that compressional wave reradiation occurs at 
nearly the same time as diffraction from the corner, so a 
Superimposed signal is received which is larger than that 
predicted by rigid impenetrable diffraction theory alone. 

For angles corresponding to backscatter from the 
leading edge of an ice lead (6 = 6, = 105°, 120°, 135°), some 
of the incident acoustic energy is again mode-converted into 
vibrations of the plexiglass wedge. As above, from an energy 
conservation point of view, less energy is available to be 
diffracted since mode conversion has occurred. However, in 
this case, compressional and shear wavesS propagate along the 
plexiglass sheets away from the corner, and hence there is no 
opportunity for reradiation from a discontinuity to be 
superimposed with the diffraction. Consequently, the measured 
diffraction is less than predicted. 

The finite size of the source and receiver may be one 
possible explanation for the excess diffraction measured at 
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6 = 6, = 120°, 135°. A true point source/receiver would be 
able to measure the actual diffracted pressure in the trough 
of the theoretical curve, whereas the finite LC-10 source/ 
receiver mounting was 2.3 cm wide at a range of 25 cm, so 
the acoustic signal was measured over a 5.3° sector centered 


at @ = 6,4 = 135°, and not at a point. 


B. BACKSCATTER FROM THE ARCTIC LEAD 

The experiments performed on the underwater 270° acrylic 
wedge (Chapter VI, Part A) and measurements made by Denny 
and Johnson (1986), indicated that during lagging edge 
backscatter, more energy was received from the edge of an 
acrylic wedge than was predicted by the Biot-Tolstoy 
diffraction theory. The difference is due to the mode 
converted energy which reradiates from discontinuities into 
the water (in some cases at the same time as diffraction 
from the plate edge). 

The objective of this experiment was to obtain 
backscatter measurements from the edge of the floating 
acrylic plate which models a smooth Arctic plate for a point 
source and receiver positioned monostatically beneath the 
plate. It iS a repeat of one experiment performed by Denny 
and Johnson (1986), but by judicious choice of source/ 
receiver positions, taking into consideration the theorerm 
cal paths and speeds of propacation for the different types 


of waves, it was possible to separate in time the received 
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signals corresponding to the diffraction, compressional, and 
shear waves which radiated from the plate edge back into the 
water for certain geometries. At the remaining angles, 
complex phasor subtraction provided the separation. 
Pen onper imental Procedure 

The laboratory measurements of the plate edge 
diffraction, and of the plate compressional and shear waves 
which reradiate from the plate edge were made in the Ocean 
Acoustics Wave Facility tank described in Chapter II. The 
source and receiver were mounted together on nylon strings 
threaded through a series of eye screws forming a pulley 
arrangement. One end of the pulley was fixed at the surface 
of the water near the center of the tank, while the other 
end was attached to the bottom of a pole positioned in a 
vertical guide along one side of the tank. Source/receiver 
positioning was then accomplished by setting the depth of 
the pole at one end, and adjusting the pulley mechanism to 
move the hydrophones along line AB, as shown in Figure 6.5. 

The floating acrylic plate was positioned above the 
source/receiver without disturbing the pulley arrangement, 
and was held ina fixed position by clamping its outer edges 
to the sides of the tank. The source/receiver was posi- 
tioned 60 cm from the plate edge, and subsequent positions 
were chosen along the same arc of 60 cm radius. The exact 
position relative to the diffracting edge was determined 
acoustically by measuring the travel time of the vertically 
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Figure 6.5: Experimental set-up for the Arctic open lead. 
Source/Receiver positioned for monostatic 
backscatter, (r = ro), (68 = 84). Two positions 
are shown. 
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propagating reflection from the plate and the travel time of 
the signal which diffracted from the edge of the plate. 
Once positioned, the source and receiver remained fixed 
while data were collected for all frequencies of interest. 

Since diffraction occurred not only from the plate 
edge of interest shown in Figure 6.5, but also from the two 
Sides and the back edge of the floating plate, it was 
necessary to subtract the effects of the three unwanted 
edges which interfered with the diffraction of interest. By 
clamping a second, identical plate to the edge of interest, 
the measurement for the same source/receiver position was 
repeated without disturbing the first plate. The received 
Signal then obtained contained no diffraction from the edge 
of interest, and diffraction from the three interfering 
edges was identical to the initial measurement. By a point- 
by-point computer subtraction of the two signals (one from 
the single plate, and one from the plates clamped together), 
the effect of the three interfering edges was eliminated, 
and only the diffraction and reradiation from the edge of 
interest remained. The quality of the computer subtraction 
depended on two crucial points: 

(a) Tf the original plate was disturbed during the 
clamping process, the position of the three inter- 
fering edges would be slightly altered, and an 
imperfect subtraction would result. 

(b) If the two plates were not clamped with the same 
depth of immersion, a discontinuity would be present 
at that joint, resulting in a small amount of 


diffraction and reradiation from the discontinuity. 
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In this case, subtraction would not provide an 
accurate indication of the diffraction from the edge 
of the original plate. 
Figure 6.6 shows a sample waveform for position Rl and the 
effectiveness of subtraction for removing reflection and 
unwanted edge interference from the signals of interest. 

The frequencies chosen for this experiment were 
those which provided a measurable diffraction from the plate 
edge, and allowed for simplified signal processing: 3 cycles 
of a 46.9 kHz sinusoid, 4 cycles of 62.5 kHz, and 5 cycles 
of 78.1 kHz. The received signal was amplified 100 times 
and filtered from 3 to 100 kHz by the Ithaco Low Noise 
Preamplifier, amplified an additional 20 dB, and sent to the 
analog to digital converter and signal processor. The 
received time-domain signal was averaged 500 times to 
eliminate random background noise, and subtracted, as 
described above, to eliminate the coherent "noise" due to 
the other three edges. The transient signals of interest 
were then transformed to the frequency domain using a 64 
usec time window and a 32-point FFT to obtain their 
magnitude, since the sampling frequency was 500 kHz. 

2. Results and Discussion 

Figure 6.7 and Table IV show the geometry and the 
individual legs of the propagation paths for six source/ 
receiver positions. Using the appropriate values for the 
speed of sound in water (cy = 1481 m/sec), compressional 
wave speed in the plate (Cp = 2257 m/sec), and the shear 
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Figure 6.6: Sample waveform for the Arctic lead experiment 
Ome = 46.9 KHz. Top picture shows the signal 
return from the single plate: specular reflection to 
the left of the cursor, edge diffraction plus 
interference to the right. Bottom picture shows the 
subtracted signal after clamping the second plate: 
specular reflection is eliminated by subtraction, 
interference is eliminated from the three unwanted 
edges leaving reradiated signal from the edge of 
interest. 


compressional wave reradiation (t = 680 usec) 
plate edge diffraction (t = 740 usec) 


(1) 
(2) 
(3) plate shear wave reradiation (t = 830 usec) 
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Geometry for the Arctic open lead experiment. 
Source/receiver monostatic backscatter posi- 
tions Rl - R6é. 


6. = critical angle for generation of plate 
compressional waves 


--- = path for optimum reradiation of plate 
compressional wave 


path for optimum reradiation of plate 
shear wave 
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SOURCE/RECEIVER POSITIONS FOR 


TABLE IV 


THE ARCTIC OPEN LEAD EXPERIMENT 


POSITION 
Rl R2 R3 R4 R5 R6 

tr (usec) 366 460 540 628 665 693 
d; (m) Oi 7 | 0.341 0.400 0.465 0.493 0151/3 
tnor (usec) = 838 836 834 = = 
te (usec) 842 oi = = = = 
tq (usec) Oat. = = = B25 821 
Ry (m) 05675 0.621 0.619 0.618 Sol Crowe 
rj; (m) 0.618 0.519 0.472 0.407 02361 02.326 
8 (deg) Be 33.3 40.3 48.8 53.8 57.5 
hy (m) Oi35.3.5 G2 422 0.494 OSS 0. 609 0.634 
1; (m) 0.421 Oy 7 1. Oe 1esa 0.069 = = 

NOLES $ 

(1) A “total" signal is present at positions R2, R3, R4 


(2) 


(3) 


where the compressional reradiation and diffraction 
overlap. 


By Snell's Law, compressional signals in the plate 
can only be generated at positions Rl - R4. 


The diffraction from the plate edge can be measured 
at position Rl where there is separation from 
compression, and at R5, R6 where there is no 
compressional interference. 
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wave speed (cg = 1373 m/sec), the predicted arrival times 
for the various propagation paths were calculated using 


Equations 6.10, 6.11) (wands6. ie. 


(Diffracted path) tp = 2Ri/cw (6.10) 
(Optimum path tc = [(hy + Rp7 cee co (6 ode 
for compression) 

(Optimum path ts = [(dy + Ri)/ceg] + [ri/cgs] (63) 


for shear) 


A sample waveform with the different separated 
arrivals is shown in Figure 6.8 where the signals due to the 
Various propagation paths are identified. 

For receiver positions R2, R3, and R4 (shown in 
Figure 6.7), the diffracted energy and the energy reradiated 
from the plate edge due to the compressional wave were 
received at nearly the same time, since the higher speed of 
propagation for the compressional wave offset its greater 
path length. By expressing each path length in terms of the 
number of wavelengths travelled by the pulse sound, the path 
difference corresponded to the relative phase shift between 
the overlapping arrivals (one wavelength is equivalent to 


360° of phase shift). Consequently, the magnitude of the 
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Figure 6.8: Sample waveform at 78 kHz for the Arctic open 
lead experiment (corresponding to position Rl 
of Figure 6.7) Starting times are: 


t = 262 usec ==> Compressional wave reradiation 
from plate edge 


300 ysec ==> diffraction from plate edge 


Gh 
ll 


t = 366 usec ==> reradiation from plate edge due 
to plate shear wave. 


(Waveform at the right-hand side of the figure 
represents interfering noise from a different path.) 
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compressional wave reradiation could be inferred from a 
vector addition of signals as shown in Figure 6.9. The 
phase difference between diffraction and compression was 
calculated from the difference in path lengths, the magni- 
tude of the diffracted signal was approximated by uSing the 
Biot-Tolstoy theory for a single rigid wedge, and the 
magnitude of the total signal was obtained by transforming 
the total received transient signal to the frequency domain. 
The results of these calculations are plotted on the polar 
diagrams of Figures 6.10 , 6.11, and 6.12. The notes 
associated with Table IV explain the presence or absence of 
experimental data points for each source/receiver position. 
For comparison, the theoretical curve plotted on each figure 
represents the Biot-Tolstoy predicted diffraction from a 
rigid wedge. 

Observation of the geometries shows that for shallow 
grazing angles for trailing edge backscatter, it is the 
compressional path which becomes superimposed with the 
diffraction; the shear path interferes with diffraction for 
trailing edge backscatter at nearly normal incidence. 
Figures 6.10, 6.11, and 6.12 clearly indicate that the Biot- 
Tolstoy rigid wedge theory, which yields the maximum energy 
for diffraction, does not account for all of the energy 
backscattered from the trailing edge of an acoustically non- 


rigid wedge. 
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Brgure 6.9: 





Wector -SuUDEract monmeorma1fttraction from the 
total signal backscattered from the edge of the 
floating plate. 


T = Total (measured) backscattered signal from 
the plate edge when compression and 
diffraction overlap (Positions R2, R3, R4) 


C = Compressional signal reradiated from plate 
edge 
D = Theorelticaieeier-ltolstey Diffraction from 


the plate edge 


¢ = Phase difference between C and D due to 
path length difference. 
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BSS due to compressional and shear waves 
reradiated from an acrylic plate floating on 
water compared to Biot-TOIlstoy predivered 
diffraction from a single rigid 270° wedge. 
Frequency, f = 46.9 kHz. Diffraction in dB ref 
spherical spreading over (r = ro = 0.6 m). 


O = TOTAL (compression + diffraction) 

0 measured compressional wave reradiation 

* = calculated compressional wave reradiation 
(vector subtraction) 

measured shear wave reradiation 

measured diffraction 

Biot-Tolstoy Theory 
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BSS due to compressional and shear waves 
reradiated from an acrylic plate floating on 
water compared to Biot-Tolstoy predicted 
diffraction from a single rigid 270° wedge. 
Frequency, £ = 62.5 kHz. Diffraction in dB ref 


spherical spreading over (r = rg = 0.6 m). 

{1} = TOTAL (compression + diffraction) 

0 = measured compressional wave reradiation 

* = calculated compressional wave reradiation 
(vector subtraction) 

A = measured shear wave reradiation 

x measured diffraction 


Biot-Tolstoy Theory 
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Figure 6 ue. 
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BSS due to compressional and shear waves 
reradiated from an acrylic plate floating on 
water compared to Biot-Tolstoy predicted 
diffraction from a single rigid 270° wedge. 
Frequency, f = 78.1 kHz. Diffraction in dB ref 


spherical spreading over (r = ro = 0.6 m). 

0 = TOTAL (compression + diffraction) 

O = measured compressional wave reradiation 

e = calculated compressional wave reradiation 
(vector subtraction) 

A = measured shear wave reradiation 

xX = measured diffraction 


Biot-Tolstoy Theory 
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C. MODE CONVERSION STUDY OF A FLOATING PLATE 
After the experiment described in Section B successfully 
generated compressional and shear waves in the floating 
plate which reradiated from the plate edge separately from 
the diffraction, the following experiment was devised to 
study the effects of plate compressional wave reradiation 
alone. Recall that, for some source/receiver positions in 
the previous experiment, the compressional reradiation and 
diffraction overlapped, and the magnitudes of each were 
identified by the vector subtraction of the arrivals. 
1. Experimental Procedure 
Figure 6.13 shows the experimental set-up. The 
source used was the PVDF plane wave transducer described in 
Appendix A, suspended beneath the floating plate at 
different angles of incidence with the acrylic surface. The 
receiver was an LC-10 hydrophone positioned beneath the 
opposite edge of the plate by an extending arm from the side 
of the tank. Diffraction from the receiver end of the plate 
was prevented by inserting a sound reflecting barrier in the 
direct path between the source and the receiver. The 
barrier was a 1" thick sheet of Celotex™, covered by a thin 
sheet of plastic to prevent water penetration and weighted 
down with 50 pounds of lead to achieve slight negative 
buoyancy. Celotex, a semi-rigid polyisocyanurate foam-board 
insulation, has closed pores filled with air, and provided a 
vertical "pressure-release" boundary beneath the center of 
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Figure 6.13: Experimental set-up for the mode conversion 
study. 


Source = PVDF planar transducer driven with all 
12 elements driven with equal amplitude and 
phase. 


Receiver = LC-10 hydrophone at positions i = 1 
- 14 (distance from plate edge = r ; depth = d) 


(NOTE: predicted paths for generation and reradiation of 


both plate compresSional waves and head waves are shown, at 
the receiver positions indicated in Table V.) 
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the plate. With this barrier in position, the direct water 
path from the source to the receiver was attenuated in 
mecess Of 50 dB, and the interfering direct and direct/ 
diffracted signals were effectively eliminated. 

The source was driven with equal amplitude elements 
to produce an early uniform plane wave. For two different 
orientations of the source, the acrylic plate was excited at 
each of the four (4) frequencies of interest, and reradia- 
tion from the plate on the opposite side of the barrier was 
measured at 14 different receiver positions, varied in both 
range and depth from the plate edge. The frequencies chosen 
Pesos >.6, 23.4, 31.3, and 39.1 kHz, using a burst of 2, 3, 
4, or 5 cycles, respectively. They were chosen both for the 
fact that the signal processing was simplified and for the 
efficiency of the plate compressional wave excitation at 
those frequencies. The angle of incidence 6; = 37° was 
chosen as one source orientation because that is the 
critical angle for generation of compressional waves in the 
aerylic plate. The angle of incidence 6; = 61° was chosen 
aS a compromise between a shallow grazing angle to model 
convergence zone propagation in the Arctic Ocean and the 
finite length of the acrylic plate which precluded very 
shallow grazing angles. 

2. Results and Discussion 

Figure 6.13 shows a simplified schematic diagram of 

the source/plate/receiver orientation with the dotted lines 
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representing the various path lengths of interest. Accurate 
knowledge of the path lengths was essential to predict the 
received signal arrival time. 

To find the effective Starting point Ofmeene 
compressional wave in the plate, the pressure profile was 
measured above the source at a depth of 1/4" below the 
plate. The 19 cm high source produced a narrow beam which 
interacted along the bottom of the plate over a 25 to 35 cm 
length, depending on the source orientation. Figure 6.14 
shows the measured pressure profile for the two source 
orientations. Since the pressure peaked at a specific 
location beneath the plate, shown on Figure 6.14, that 
distance from the plate edge was assumed to be the point of 
maximum excitation of the plate compressional wave, shown as 
point B on Figure 6.13. 

At the receiver end of the plate, it was predicted 
that the presence of the plate compressional wave would be 
observed either by reradiation from the edge of the plate as 
in the previous experiment, or as head wave radiation. Head 
waves are described by Clay and Medwin (1977) and are well 
known by geophysicists for a thick, high-speed bottom 
underlaying the lower speed ocean water. A similar situa- 
tion appears to exist for the high speed compressional wave 
in the Arctic ice overlaying the lower speed ocean water. 
Head waves are generated when a sound signal travelling in 
one medium refracts into a second, higher speed medium by 


82 


Source pressure 1¢m below plete (d& ref max) 


0 


—-15 





Figure 6.14: 


Source pressure vs position 


15 20 
Distance along plate, x 


source pressure vs. 


3 


es 


(cm) 


position. 








the critical ray path, causing bulk waves to propagate along 
the interface at the higher sound speed. The presence of 
these "lateral waves" is indicated by their reradiation at 
the critical angle back into the lower speed medium. 

Table V lists the predicted arrival times for the 
head wave and for compressional wave reradiation at the 
plate edge, as shown in Figure 6.13. The measured arrival 
times agreed with the path for plate edge reradiation for 
receiver positions 2, 7, 13, and 14, since those positions 
were outside the cone of head wave reradiation. All other 
arrival times agreed more closely with the path for head 
wave reradiation, and the plate edge radiation was masked by 
the stronger head wave. Consequently, it was impossible to 
obtain the desired radiation pattern from the plate edge for 
which this experiment was designed. However, valuable 
information was obtained by discovering the presence of the 
head wave. 

Figure 6.15 shows a plot of the head wave spectral 
density amplitudes in the relative receiver positions for 
each of the three source frequencies with the source 
oriented at 6; = 61°. Although the measured amplitude of 
the head waves do not precisely agree with the theoretical 
behavior described in Clay and Medwin (1977) given as 

amplitude of head wave = ace ety 


the character of the reradiation is correct since the 
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TABLE V 


PREDICTED AND MEASURED ARRIVAL TIMES 
FOR PLATE COMPRESSIONAL WAVE RERADIATION 


Source Angle 6; = 37° 


Arrival time (usec) 








Receiver Predicted Predicted Actual 
Position d(m) r(m) Head Wave Plate edge 
1 0.05 0.02 716 728 710 
Z 0305 0 20 762 767 761 
2 Dose lind: 5 1, 707 792 675 
4 Cert -0.05 ae TT 105 
5 Ovaikel: 0.00 72s 766 729 
6 On! 0.05 7959 TiS yo 
i Ons 0.10 oe) 792 780 
8 0.16 -0.20 709 864 664 
9 O26 -0.14 wey, 835 690 
10 0.16 -0.07 22 810 718 
il. 0.16 O20 0 743 800 T37 
12 0.16 0.07 786 810 800 
ee 0.16 0.14 830 835 826 
14 0.16 0.20 859 864 867 
Source Angle 6; = 61° 
Arrival time (usec) 
Receiver Predicted Predicted Actual 
Position d(m) r(m) Head Wave Plate edge 
1 0.05 0.02 74] 767 719 
2 0.05 0.10 794 809 785 
S Oe rt -~0.10 vou 837 691 
4 OL -0.05 WS 821 720 
5 On): 0.00 759 810 743 
6 Ou ch 0.05 796 821 712 
7 Omi 05.50 811 837 796 
8 0.16 -0.20 roo 909 680 
3 EG -0.14 oa 880 704 
10 0.16 -0.07 756 854 741 
ia 02256 0.00 TAL 844 768 
i 2 0.16 0207 818 854 802 
ie 0.16 0.14 864 880 844 
14 O26 O20 894 909 879 


(NOTE: Actual arrival times represent the average arrival 
time at that receiver position for the four frequencies of 
interest.) 
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amplitudes approximately follow lines of constant pressure. 
For the parameters given above, x is the horizontal distance 
between source and receiver, L is the distance that the head 
wave travels in the plate, and k is the wave number. 

The source waS oriented to correspond with the 
critical angle ray path (8; = 37°), as well as a shallower 
angle (6; = 61°) at which total internal reflection would 
also take place, to simulate how sound propagating over long 
ranges in the Arctic will interact with the ice cover (such 
as for convergence zones). Mode conversion of the compres- 
sional wave in water to the compressional (or lateral) wave 
in the plate was 40% more efficient when the source was 
oriented at the critical angle, as indicated by the head 
wave magnitudes for 15.6 kHz in Table VI. The behavior of 
the head waves was similar for the other frequencies, 
Suggesting that for convergence zone propagation in the 
Arctic, both mode conversion and reflection may occur at 
each interaction with the undersurface of the ice, since 
head waves were also effectively generated at shallow 


grazing angles. 
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Figure 6.15: Head wave amplitudes (rounded to nearest 
integer) for mode-conversion experiment at each 
frequency. Circled values at bottom of each 
figure represent approximate location of 
isopressure lines. Receiver positions marked by 
an 'x' indicate those positions outside of the 
cone of head wave radiation. 


87 


TABLE VI 


MEASURED HEAD WAVE AMPLITUDES 
(for frequency, f = 15.6 kHz) 


Receiver Source Angle 6; = 37° Source Angle 6; = 61° 
Position Amplitude, (V_ pp) _ Amplitude, (V_pp) 
il 2 Gus OU E>. 35 
2 9.99 7.49 
3 27.00 ee? 
4 32250 2067 
5 18.43 PE. 02 
6 23.68 1:4), 055 
y S356 MOKG5 
8 40.54 Zire US 
9 40.91 19.42 
10 34.19 19.14 
ie 22.48 Foe 3 
2 22 14.07 
13 19.93 oR 2 
14 13.34 hei 
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Wid ihe ELEVATION ANGLE DEPENDENCE OF 
BACKSCATTER FROM A RIDGE 
The objective of this experiment was to obtain backscat- 
ter measurements from the three acrylic models of different 
size ice keels in order to determine how the level and the 
character of backscatter depend on the height and slope of 
the keel. Measurements were made to show the characteristic 
behavior of backscatter for grazing and normal incidence, as 


well as the dependence on frequency and elevation angle. 


Pee SAPERIMENTAL PROCEDURE 

The laboratory measurements of backscatter from the ice 
pressure ridge models were performed in the Ocean Acoustics 
Wave Facility tank described in Chapter II. The fabrication 
process and dimensions of the three pressure ridges are 
described in Chapter III. 

Two LC-10 hydrophones were mounted for monostatic 
backscatter by extending a single, thin stainless steel rod 
from the side of the tank. Typical measurements were made 
with the source/receiver positioned 50 cm from the "ice 
keel", which prevented plate side diffraction from interfer- 
ing with the backscattered signal of interest, and also 
allowed the strong direct signal between source and receiver 


to decay prior to receiving the backscattered signal. 
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At shallow grazing angles, the reftlecticn i romeene 
overhead plate was well isolated and could be windowed out 
from the ridge scatter. At steep grazing angles (greater 
than 50°), the surface reflected signal from the smooth 
plate directly above the source/receiver, overlapped the 
backscattered signal of interest. Figure 7.1 illustrates 
the technique used to overcome this difficulty. Figure 7.1 
(a) is a typical waveform for an incident angle of 50°, and 
shows the strong surface reflection immediately followed by 
the signal backscattered from the ridge. Moving the plate 
to increase the propagation path to the ridge results in 
Figure 7.1 (b), where the backscattered signal has been 
shifted out of the time window and only overhead reflection 
remains. Figure 7.1 (c) is an ineffective point-by-point 
computer subtraction of Figures 7.1 (a) and 7.1 (b); 
however, to obtain backscatter only, a slightly different 
depth of plate immersion resulted in a measurable time shift 
in the reflected arrival sufficient to prevent clean 
subtraction. Figure 7.1 (d) is the resulting subtracted 
waveform when the time shift was corrected prior to subtrac- 
tion, and only the backscattered signal remains. 

Four experiments were performed on each of the three 
pressure ridges: 

(1) A horizontal probe (at the depth of the ridge crest) 


of the backscattered pressure field due to the ridge 
at grazing over the range r = ro = 20 to 70 cm. 
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Subtraction technique used in backscat- 
ter measurements. Grazing angle = 50°, 
frequency = 62 kHz. 
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Plate reflection (1) plus total 
backscatter (2). 

Plate reflection alone. 
SDM Of Cajee— (b) without 
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variation in plate immersion depth. 
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usec time shift to account for the 
difference in plate immersion 
depth. The major signal is due to 
ridge backscatter 
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(2) A vertical probe of the backscattered pressure field 
below the ridge for normal incidence over a depth 
range d = dg = 10 to 70 cm. 

(3) A frequency dependence experiment for which the back- 
scatter was measured for nine frequencies at grazing 
angles 6 = 10° and 6 = 20° from the horizontal: 

(4) An experiment to determine the dependence of back- 
scatter on elevation angle, by holding the source/ 
receiver at a fixed radius r = ro = 508cm an@ prema 
in elevation angle from 6 = 0° to 6,4 = 90° to find 
the "directivity pattern” of the sound scattered from 
the ridge. 

Figure 7.2 shows the source/receiver geometry with 
respect to the ridge for each of the four experiments. 
Typical measurements were made driving the source with a 20 
volt peak-to-peak pulsed sinusoid of l, 2 or 3 cycles and 32 
usec duration. Although the electrical input signal was 32 
usec, the received acoustic signal was much longer due to 
mode-converted reverberation in the ridge, and backscatter 
from different parts of the ridge. 

The analysis window chosen was sufficiently long to 
include all of the backscattered energy, then "zero-padded" 
with as few extra zero points as possible in order to obtain 
a power of two (2 total points, n = any integer) to allow 
the received signal to be transformed uSing a Discrete 
Fourier Transform (DFT) algorithm. The source frequencies 
used coincided with the bin frequencies from the DFT to 
minimize the Gibbs' effect (Clay and Medwin, 1977). The 
backscattered spectral pressure from the ridge was propor- 


tional to the magnitude of the response from the DFT, and 
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all received signals are expressed in decibels, using the 


backscattering strength (BSS) definition: 
2 -_N + 2 
Prac ae | 
Pawnee ! 
aa 


In this equation, Pg is the reference pressure at range Ro, 


al 


BSS (dB) = 10log,, 


PeS is the received backscatter pressure at range R, A is 
the surface area over which the scattering 1S measured, and 


N is the degree of geometrical divergence of the backscatter 


The Biot-Tolstoy theory of diffraction predictema 
maximum amount of energy which will diffract from an 
infinite, rigid wedge. When digital signal processing 
technigques are used, the wedge need not be infinite in 
length if the analysis time window is closed before sound 
reaches the ends of the finite wedge. Naturally, if the 
wedge is not infinite in length, something less than the 
maximum amount of energy will be diffracted. The cal- 
culation of data for Figure 7.3 was provided by Mr. Ken 
Reitzel and shows the amount of energy which will diffract 
from various length finite wedges in comparison to the 
infinite wedge, where the wedge length is 2Y, centered about 
the least time point. Notice that only a fraction of the 
maximum possible energy is diffracted when 2Y is less than 


an acoustic wavelength, and nearly all of the energy 1s 
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diffracted when 2Y is greater than a wavelength. A finite 
wedge whose length equals one acoustic wavelength will 
Giffract 96% of the maximum possible energy. From another 
point of view, 96% of the diffracted energy comes from one 
acoustic wavelength of the ridge in the least time region, 
and the remainder of the ridge contributes only 4% of the 
energy. For backscatter at long ranges, a pressure ridge 
should effectively behave as a point scatterer (since the 
majority of the energy is diffracted from the least time 
region), and spherical divergence would be expected. Then, 


the equation for backscattering strength becomes: 
2 4 
2 acl a 
BSs (dB)y = LORegrs aoe i (7 ae 
Je 


Ir. the case of the vertical probe experiments, the 
backscattered signal was essentially a reflection from the 
smooth plate, slightly altered by the presence of the 
ridges. In that case, the backscatter can also be calcu- 


lated in terms of dB referenced to a mirror (1.e., perfect 


reflection): 


| Pag R 
GB (ref mirror) = Z0logay re (7.28 
Po Ro 
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If the backscatter comes from a point, then experimental 
agreement with theory will be better when Equation (7.1) is 
used; if the backscatter more closely resembles a reflec- 


tion, then Equation (7.2) yields the best agreement with 


theory. 


5. "TYPICAL" RIDGE RESULTS 

The "typical" ridge used for this experiment was 
designed to model the typical pressure ridge found in the 
Central Arctic basin as described in Chapter III. It served 
as the standard against which the "large" ridge (twice the 
height) and the "steep" ridge (twice the slope) were 
compared. 

1. Horizontal Probe 

For the horizontal probe using the "typical" ridge, 

the source/receiver were positioned 1.6 cm below the ice 
plate to correspond with the depth of the ridge crest. 
Figure 7.4 shows a typical received waveform. The received 
Signals in all cases were composed of the two parts shown: 
an initial arrival which agreed in time with diffraction 
from the ridge crest (corresponding to point B) and a second 
arrival which agreed with the propagation time required for 
an acoustic signal to be mode-converted into the ridge, 
reflect from the back facet of the ridge (at point C), and 
return by the reverse path. There was no measurable 
backscatter observed from point A of the ridge. 
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PIigune 7.4: 
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Typical received waveform for horizontal 
probe using the "typical" ridge. 


(f = 62.5 kHz, N = 2 cycles, r = 40 cn, 
500 mV scale) 


(1) = diffraction from the ridge crest at 
Pelrnmees 


(2) 


reflection of mode-converted energy 
From DPoume Cc 
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Figure 7.5 shows the divergence characteristic for 
Maes nOrizZontal probe from the "typical" ridge. Equation 
(7.1) was used to calculate the BSS, assuming the 
backscatter from the ridge diverged spherically. Except for 
oscillations, the BSS is virtually constant with increasing 
range, suggesting the ridge must have produced the backscat- 
ter from an effective point (corresponding to the least time 
region on the ridge) for spherical divergence. Since the 
source/receiver were quite close to the plate, the 
oscillations of BSS with increasing range are probably due 
to surface interference of the backscattered sound. 

2. Vertical Probe 

In the vertical probe experiment on the "typical" 
ridge, certain effects of the characteristic behavior of the 
backscatter were apparent, but could only be analyzed 
qualitatively. It was more obvious in the vertical probe 
experiment for the "large" ridge, (where signal separation 
in time was improved due to the larger dimensions), what the 
various parts of the backscattered signal were due to. 

Figure 7.6 shows a typical received waveform for the 
case where 2 cycles of a 62 kHz sinusoid are used. The 
direct signal is slightly more than 2 cycles due to hydro- 
phone ringing. There appeared to be three distinct com- 
ponents to each backscattered waveform. The first arrival 


in time agreed with the path length for diffraction from the 
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Figure 7.5: Plot of BSS versus range, r (cm), for the 
horizontal probe experiment using the "typical" 
ridge. (mean value of BSS in parentheses) 
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Typical received waveform for vertical 


probe experiment using the "typical" 
ridge. 


(f= G2 oe pe Nee eC lLeceerd = 40cm, 
2.0 V_ scale) 


(l) = diffractvenerrem- crest of ridge 
(path S-A-R) 


(2) ="seflect won. tf romerce-alireinterface 
(path S-A-B-A-R) 


(3) = reradiation of ridge reverberation 
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crest of the ridge. The second; snuch larger arrival agreed 
in time with the path length for an acoustic. signal 
(13 usec later) which reflected from the "ice"-air interface 
before returning by the reverse path. Finally, theme 
appeared to be some energy arriving after the reflection, 
due to the reverberation of mode-converted sound within the 
ridge which was then reradiated into the water. 

Figure 7.7 is a plot of the BSS versus depth below 
the ridge for the vertical probe, calculated using Equation 
(7.1). Since the 31 kHz BSS was nearly independent of depth 
for range greater than 35 cm, the assumption of spherical 
divergence from a point on the ridge is evidently valid. At 
62 and 93 kHz, the variations of the backscatter about the 
average value is assumed to be due to interference of the 
spatially-extended ridge reverberation and reflection from 
the acrylic-air interface. When the scattering is plotted 
in terms of dB (ref mirror) as in Figure 7.8, the results 
were also nearly linear, which implies that the effects of 
the plate reflection were just as important as the backscat- 
ter from the ridge since the experimental agreement with 
theory is about the same when using either definition. 

3. Frequency Dependence 
The frequency dependence experiment using the 


"typical" ridge produced no meaningful results. 
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megane 7.7: Weakot ofw BSS versus depth, d(cm), for the 
vertical probe experiment using the "typical" 
ridge. (mean value of BSS in parentheses) 
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4. Elevation pice Dependence 


The source and receiver were positioned for monos- 
maere backscatter at a fixed distance r = ro = 50 cm from 
the ridge for grazing angles ranging from 0° to 90°. #£=Two 
source frequencies were used, f = 62.5 kHz and f = 93.9 kHz, 
corresponding to A, = 2.37 cm and ee Oeil. 

Assuming that the ridge reverberation would radiate 
from the facets, the ridge was modelled as a finite reflec- 
tor or aperture, Similar to a rectangular piston, infinite 
maeone direction. In the case of the "typical" ridge, the 
"active" piston face, corresponding to the slant face 
height, was H = 3.41 cm so that 

H/, = 1.44 BN. = 2.16% 
The directivity pattern for the equivalent rectangular 


piston transducer is given by: 


k 
sin(“3) sine ates 


hae 


— 


eo, k) = a 7 
(Ss 


) sind ( 


ss) Sine 

where ¢ and y are the angles shown in Figure 7.9, and W and 
H are the width and height of the piston. All measurements 
were taken for source/receiver positions normal to the 
Mmagiae, so that ¢ = 0° and 0° $ y $s 90°. The magnitude of 
the backscattering strength was plotted as a function of 


angle and compared to the theoretical directivity pattern 


which was adjusted to a level for best fit as shown in 
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Pigunes 7.20 Geometry for radiation from a rectangular 
piston. (Clay and Medwin, 1977) 
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Figure 7.10. The piston axis, xy = 0, was set to be perpen- 
dicular to the facet. 

The backscattered data were straightforward to 
analyze. Backscatter measurements for grazing angles 6 = 0° 
to 50° could be analyzed directly since there were no 
interfering signals. For grazing angles, § = 60°, 70°, and 
80°, the specular reflection from the overhead plate 
interfered with the backscattered signal and was subtracted 
out using the technique described in the experimental 
procedure section. At normal incidence (9 = 90°), the plate 
reflection and backscatter completely overlapped (as in the 
vertical probe experiment), resulting in a measured value 
much larger than predicted. It would be meaningless to 
subtract the reflected portion of the backscattered signal 
at normal incidence because it is an important part of the 
scattering characteristic for that geometry, as shown in 
mroure 7.10. 

There was fair agreement between the theoretical 
curve and the measured backscattering strength as a function 
of angle for the "typical" ridge. The BSS measurement has 
two components: the Biot-Tolstoy diffraction from the three 
corners of the ridge, as well as the radiation from the 
rectangular apertures. Since the characteristic dimension 
of the piston was on the same order of magnitude as an 


acoustic wavelength, the rectangular aperture was not an 
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efficient projector, so neither effect dominated. Cence. 
quently, the comparison with the equivalent piston radiation 
pattern was incomplete, since diffraction from the corners 


also represented a Significant contribution to the BSS. 


oF "LARGE" RIDGE RESULTS 

The "large" ridge was twice the height but the same 
slope as the "typical" ridge (see Figure 3.1 and Table III 
for the dimensions). The experiments using the "large" 
ridge provided keen insight into the response of an ice keel 
to an incident acoustic wave; since the dimensions provided 
greater separation in time of the signals backscattered from 
various parts of the ridge, and therefore, allowed each 
effect to be analyzed individually. 

1. Horizontal Probe 

For the "large" ridge, the source/receiver were 
positioned 3.2 cm beneath the plate to correspond with the 
depth of the ridge crest. 

Figure 7.11 shows a typical received waveform and 
the probable paths which combined to yield those results. 
Notice that diffraction from the crest of the ridge repre- 
sented the largest contribution to the backscattered signal. 
When the incident acoustic energy struck the front facet of 
the ridge, a small amount of energy was diffracted from the 
Sesconcinuity at point A. Then, the incident wave reached 
the ridge crest (point B) and a larger amount of energy was 
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diffracted. Finally, some of the energy incident on the 
front facet was mode-converted into compressional waves in 
the ridge, propagated through and reflected from the back 
side of the ridge (point C), and returned to the receiver by 
the same path to complete the backscatter. Bach of the 
paths described above has been verified by the arrival time 
using the specific path lengths and appropriate speeds of 
propagation. 

Figure 7.12 shows a typical divergence characteris- 
mre fOr a horizontal probe from the ridge. For spherical 
divergence of sound from the point source, and spherically 
diverging backscatter from the ridge, overall divergence of 
pressure proportional to the inverse square of the range was 
expected. BSS was calculated using Equation (7.1), where 
spherical divergence was assumed. The experimental results 
were in excellent agreement with the inverse square be- 
Gemror, indicating that the majority of the backscatter 
occurred from the region on the ridge near the least time, 
and that the backscattered divergence was spherical. 

2. Vertical Probe 

In the vertical probe experiment on the "large" 
ridge, the size again provided relatively long time separ- 
ations between the different propagation paths to allow each 
portion of the total backscattered signal to be analyzed 


iadividually. 


Lash 
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FIGULe 7 ce Plot of BSS versus range, r (cm), for the 


horizontal probe experiment using the 


"large" ridge. (mean value of BSS shown 
in parentheses) 
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Figure 7.13 shows a typical received waveform and 
the probable paths which combined to yield those results. 
The total backscattered signal is the composite of three 
emstcinct Contributions: (1) diffraction from the crest of 
the ridge at point A; (2) reflection from the "ice"-air 
interface due to the energy which was mode-converted into 
the ridge, reflected from the surface, and reradiated back 
into the water (path ABA); and (3) mode-conversion induced 
reverberation in the ridge, which is seen in the waveform as 
the decaying energy after the surface reflection and 
includes the "double reflection" within the ridge (path 
ABABA) . 

Figure 7.14 is a plot of the divergence characteris- 
tic due to the diffraction from the ridge crest (correspond- 
ing to the energy arriving prior to the reflection in Figure 
7 ee Since the diffracted signal experiences no inter- 
ference from other paths, the expected spherical divergence 
was obtained, and is shown as the inverse square behavior 
for the round trip energy. 

When the total backscattered signal (diffraction plus 
reflection plus reverberation) is plotted as a function of 
position, as in Figure 7.15, the inverse square dependence 
is not obtained due to interference effects. For example, 
the large difference in the magnitude of the received 


pressure for a frequency of 93 kHz at 30 and 35 cm below the 
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Typical received waveform for fe 
vertical probe experiment using the 
"large" ridge. 


(f = 93.9 kHz NO@=e cycles, d= 30cm, 
200 mV scale) 


l1 = diffraction from ridge crest 


(path S-A-R) 


2 = reflection from ice-air interface 


(path S-A-B-A-R) 


3 = double reflection and reradiation 


(path S-A-B-A-B-A-R) 
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eirgure 7.14; Plot of the divergence characteristic due 


to diffraction from the crest of the 
"large" ridge as a function of source/ 
Beceiver depene (The measured values 
correspond to energy received prior to 
the weet leceton win rweure 7Fak3.) (mean 
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backscattered energy.) (mean value of BSS 
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ridge is made understandable by Figure 7.17, in which the 
upper waveform is the received signal for d = dg = 30 
cm, and the lower waveform for d = dg = 35 cm. The charac- 
ter of the diffraction and the double reflection are nearly 
identical, but the reflected Signal is superimposed with the 
ridge reverberation and can result in constructive or 
destructive interference. Figure 7.16 shows the backscatter 
plotted as dB (ref mirror) for the "large" ridge, and the 
results are not Significantly different than in Figure 7.15. 
At 31 kHz, Figure 7.16 is nearly flat at long ranges 
implying that reflection is the dominant effect at low 
frequency, while at 62 and 93 kHz, the curves are not 
linear, which implies that backscatter from the ridge is 
just as important as reflection from the smooth plate at 
higher frequency. 
3. Frequency Dependence 

The frequency dependence experiment uSing the "large" 
ridge produced the only meaningful results. The backscatter 
was measured for nine frequencies at grazing angles 6 = 10° 
and 6 = 20° from the horizontal. Typical received wave- 
forms had the same characteristics as the horizontal probe 
measurements using the "large" ridge (Figure 7.18). 

Since the source and receiver were at shallow grazing 
angles and a range of 50 cm from the ridge crest, the back 


corner of the ridge (point C on Figure 7.18) was in the 
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Comparison of two waveforms for the 
vertical probe using the "large" ridge. 
The upper figure is the same as Figure 
7.13, ad = 30 cm. The lower figure (d = 
35 cm) shows how the ridge reverberation 
destructively interferes with the "“ice"- 
air reflection. 
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PIGUYreC Paneer Typical received waveform for the 


frequency dependence experiment uSing the 
"large" ridge. 
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geometrical shadow. The backscatter consisted primarily of 
diffraction from points A and B, and the plot of BSS versus 
frequency, Figure 7.19, showed an interference pattern due 
to those two backscattered components. For this source/ 
receiver geometry, the behavior of the BSS can be explained 
by accounting for the phase difference between the two 
Gentributions: 
phase difference = kH/cos6 = nT, 
where k is the wave number, H = the difference in path 
length from the source/receiver to point A or B, and 6 is 
the grazing angle. Figure 7.20 shows a plot of BSS versus 
phase difference, kH/(Tcos6), and it becomes apparent that 
the peaks occur approximately for even multiples of 7 
(constructive interference), while the troughs occur at odd 
multiples of ™ (destructive interference). 
fee elevation Angle Dependence 

The geometry, source frequencies and modelling 
technique for the "large" ridge were exactly the same as 
those described for the "typical" ridge. In this case, 
however, the height of the ridge facet, H = 6.82 cm, 


resulted in height to wavelength ratios of: 
H/, = 2.9 H/dA, = 4.32. 
The magnitude of the backscattering strength was again plot- 


ted as a function of angle and compared to the theoretical 
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Pugune= 7 eo: Plot of BSS versus frequency, f(kHz), for 
the frequency dependence experiment using 
the "large" ridge. (r= 50 cm, 6 = 10°, 
Ze) 
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Figure 7.20 (b) 


Begure 7,20: Plot of BSS versus phase difference 
kH 
n =-!|!-_ 
Incosé! 
for the frequency dependence experiment 
using the "large" ridge. 
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directivity pattern for the facet acting as a rectangular 
piston transducer, as shown in Figure 7.21. 

Since the height of the piston face was relatively 
large compared to a wavelength for the frequencies used, the 
experimental results were in excellent agreement with the 
theoretical directivity pattern. Biot-Tolstoy diffracwren 
occurred from the three corners of the ridge, but was small 
in comparison to the radiation of mode-converted reverbera- 
tion from the larger rectangular apertures. For the "large" 
ridge, note that the value of backscatter for 96 = 90° shows 
good agreement: since the ridge was large in width, it 
scattered effectively so that the measured backscatter can 
be attributed primarily to the ridge diffraction instead of 


reflection from the ice sheet behind the ridge. 


ee "STEEP" RIDGE RESULTS 
The "steep" ridge used for this experiment was designed 
to model a possible pressure ridge which could be found in 
the Central Arctic basin. It was the same height, but twice 
the slope as the typical ridge, so it was designed to show 
the dependence of backscatter on ridge slope. 
1. Horizontal Probe 
For the horizontal probe using the "steep" ridge, 
thesource/receiver were positioned 1.6 cm below the ice 
plate to correspond with the depth of the ridge crest. 
Figure 7.22 shows a typical received waveform. The received 


124 


(°wo 9G = I ‘UMOYS 2Te satouenbaly OMI) 


IOoJ stxe uoystd ayy smoys MOTTe SUL 


‘qyooey esbpta yors 
-uoqstd aetnbuej,de1 ZuaTeaAtnbs ue 


woTy uoTReTpeA YaoeReF FO uzaq3zed AATATIOSITpP 9Yy OF pareduoo abpti ,woODIeT, 


544 403 YyRbuezqys AaZQZeOSyoeq Te2IO 94 


“er 06-6 096 


9° 


7% 
way 


O A 


Ue =6 


ZHY SE 


j 





ap Ol- O2- O¢- A os 


“ ONS 


Og - 





Jo aouaepuddap aTbhue uot zeaAolTY 


[TZ°L4 eanbty 


OL =6 


on 


[ZH 29 | 


Ol - 


122 


AIR Ace aale 


A ane 


——— So 
WATER aes 
S/R ee SS 


Figures 72: Typical received waveform for the 


horizontal probe experiment using the 
"steep" ridge. 


(fo= 62.5 kH2,QNe= 22° eveles eee 940 cre 
500 mV scale) 


(1) = Backscatter from front facet (AB) 
(2) = Mode conversion and reflection from 
back facet (BC) 
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Signals in all cases were composed of the two parts shown: 
an initial arrival which agreed in time with reflection or 
diffraction from the ridge facet closest to the source/ 
receiver; and a second arrival which agreed with the 
propagation time required for an acoustic signal to be mode- 
converted into the ridge, reflect from the back facet of the 
ridge, and return by the reverse path. Due to the small 
dimensions of the "steep" ridge and the positioning error of 
the source/receiver, it was impossible to positively 
identify which part of the ridge caused each backscattered 
component, so the backscatter was generalized to the front 
facet and the back facet contributions. The signals had 
Similar characteristics to those for the horizontal probe 
from the "typical" ridge, with the difference that the 
backscatter from the near and far facets were separated by a 
smaller time interval in the case of the "steep" ridge. 
Figure 7.23 shows the divergence characteristic for 
the horizontal probe from the "steep" ridge. Equation (7.1) 
was used to calculate the BSS, assuming the backscatter from 
the ridge diverged spherically. The curve of BSS versus 
range shows oscillating and decreasing BSS with increasing 
range for the "steep" ridge. The source/receiver was the 
same depth below the plate as for the "typical" ridge 


horizontal probe, so both show similar surface interference 
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effects. The average values of BSS are lowest for the 
"Steep" ridge and highest for the "large" ridge. 
2. Vertical Probe 

In the vertical probe experiment on the "steep" 
ridge, fewer effects of the characteristic behavior of the 
backscatter were apparent than for the "typical" and the 
"large" ridge vertical probes. Figure 7.24 shows a received 
waveform. Notice that the majority of the signal appeared 
to be reflection from the ice plate. There was no 
measurable diffraction from the crest of the ridge, and 
there was significantly less reverberation from the ridge 
reradiating after the reflected signal. 

Figure 7.25 is a plot of the BSS versus depth below 
the ridge for the vertical probe, calculated uSing Equation 
(7.1). Since the BSS was nearly independent of depth, the 
assumption of spherical divergence from a point on the ridge 
must have been valid. In this case, Since reverberation was 
not significant, there was less variation of the backscat- 
tered values about the mean than was observed for the 
"typical" ridge. It is noteworthy that the plot of BSS 
versus depth decreases for shallow depths. Using Equation 
(7.1) to calculate BSS presupposed that the divergence would 
be spherical. BSS values less than the average for the 
shallow depths indicates that the divergence was over- 


compensated for, and that the actual divergence for the 
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vertical probe experiment using the 
"steep" ridge. 
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mergure 7.25: Plot of BSS versus depth, dad (cm) for the 
vertical probe experiment using the 
"steep" ridge. (mean value of BSS in 
parentheses) 
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"steep" ridge was probably dominated by reflection. When the 
scattering was calculated using Equation (7.2) and plotted 
in terms of dB (ref mirror) as in Figure 7.26, experimental 
agreement with theory was excellent: the small size of the 
"steep" ridge resulted in scattering that was dominated by 
the reflection from the smooth plate. 

3. Freguency Dependence 

The frequency dependence experiment using the 
"steep" ridge produced no meaningful results. 

4. Elevation Angle Dependence 

The geometry, source frequencies and modelling 
technique were exactly the same as those described for the 
"typical" ridge. In this case, however, the height of the 
ridge facet, H = 1.92 cm, resulted in height to wavelength 
EALYOS OF: 

H/A, = 0.81 Hee ere ee 
The magnitude of the backscattering strength was again 
plotted as a function of angle and compared to the theoreti- 
cal directivity pattern for the facet as a rectangular 
piston transducer, shown in Figure 7.27. 

The theoretical directivity pattern and the measured 
backscattering strength as a function of angle had their 
poorest agreement in this case of the "steep" ridge. The 
characteristic dimension of the piston transducer was 


approximately equal to the wavelength, so the piston behaved 
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in a highly non-directional manner. The radiation from the 
small rectangular apertures was inefficient, and Biot- 
Tolstoy diffraction which was not accounted for in the 
theoretical comparison was probably significant. 

The backscattered data were analyzed as before, with 
interfering reflections subtracted at nearly normal grazing 
angles. At normal incidence (6 = 90°), the plate reflection 
and ridge backscatter again overlapped, resulting in a 
measured value much larger than predicted. Just as in the 
vertical probe experiment for the "steep" ridge, the 
effective backscatter at normal incidence is insignificant 


compared to the plate reflection. 


ive SUMMARY OF RIDGE RESULTS 

Based on the results of the backscatter measurements for 
the four experiments on each of the three ridges, the 
following generalizations can be made. 

1. Horizontal Probe 

For all three ridges, the plot of BSS versus range 

generally obeyed a spherical divergence characteristic, 
which indicated that the majority of the backscatter 
occurred from an effective point on the ridge in the least 
time region. The "large" ridge produced the highest BSS and 
the best agreement with spherical divergence. The "typical" 
and "steep" ridges were both half the height of the "large" 
ridge, so the source/receiver was positioned closer to the 
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plate for their horizontal probes, and surface interference 
problems degraded the agreement with the spherical diver- 
gence characteristic in those cases. 
2. Vertical Probe 

Received waveforms for vertical probe experiments 
typically consisted of three components: Biot — Tome trey 
diffraction from the ridge crest; mode-conversion and 
reflection from the "“ice"-air interface; and a double 
reflection within the ridge prior to reradiation. The 
"large" ridge resulted in the greatest separation in time 
between those three components, so each could be analyzed 
individually; the "typical" ridge showed evidence of each 
component, although not well separated in time; and the 
"steep" ridge only showed signs of reflection, since the 
size of the diffracted and double reflected signals were 
negligible. The backscatter from all three ridges generally 
obeyed a spherical divergence characteristic beyond some 
critical depth below the ridge. Variation in the divergence 
characteristic at different source/receiver depths was 
explained by the destructive interference of the reradiation 
of mode-conversion induced ridge reverberation with the 
reflected portion of the backscatter, as shown in Figure 
he Mees 

With the source/receiver positioned below the ridge, 


the received signals were a combination of sound scattered 
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from the ridge, as well as the reflection from the smooth 
plate. When scattering from the ridge is most important, 
Equation (7.1) for the backscattering strength is a more 
appropriate comparison with theory; when reflection from the 
smooth plate is important, Equation (7.2) for dB referenced 
to a mirror should be used. In the case of the "large" 
ridge, scattering was just aS important as reflection, so 
experimental agreement with theory was nearly the same when 
using either Equation (7.1) or Equation (7.2). At the other 
extreme, the small size of the "steep" ridge meant that 
reflection was much more important than backscatter. As a 
result, experimental agreement with theory was better when 
plotted in terms of dB ref. mirror. 
3. Frequency Dependence 

The frequency dependence measurement of the "large" 
ridge produced the only meaningful results. For the shallow 
grazing angles of this experiment, diffraction occurred from 
two corners of the ridge which produced an interference 
effect due to their slightly different path lengths in 
comparison to a wavelength. The "typical" and "steep" 
ridges produced inconsistent results, possibly due to the 
increased surface interference problems associated with 


those two ridges for their smaller height. 


le 


4. Elevation Angle Dependence 


Each of the ridges was modelled as a finite reflec- 
tor or aperture with a height corresponding to the ridge 
slant face. The backscatter consisted of the reflection and 
reradiation of sound energy from the rectangular aperture, 
as well as diffraction from the corners of the ridge. The 
BSS was compared to the directivity pattern for an 
equivalent rectangular piston. Experimental agreement with 
theory was good for the "large" ridge because radiation from 
its large aperture (in comparison to a wavelength) dominated 
the diffraction from the ridge corners. Agreement with 
theory in the case of the "typical" or "steep" ridge was not 
as good since the diffraction increased in relative impor- 
tance as the radiation became less efficient for these 
smaller rectangular apertures. The BSS was greater for 
larger ridges: doubling the slope of the "typical" ridge 
reduced the BSS by 10 dB, while doubling the height of the 
"typical" ridge increased the BSS by 5 dB. 

Finally, it was observed that diffraction theory did 
not account for all of the energy backscattered from a 
ereigiee During the frequency dependence experiment for a 
grazing angle of 20°, there were no surface interference 
effects and no interference of specular reflection from the 
smooth plate, so the geometry was optimized to measure 


diffraction alone. The additional energy was due to the 
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reradiation of mode-conversion induced reverberation from 


the ridge, as shown in Figure 7.28. 
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VIII. CONCLUSIONS 


Due to the physical characteristics of the laboratory 
acrylic model and the ability to scale the laboratory 
dimension-to-wavelength ratio to match typical Arctic 
dimension-to-wavelength ratios, it is possible to study the 
low frequency scattering properties of the Arctic ice canopy 
under controlled laboratory conditions. These experiments 
were designed to study the behavior of two particular 
scattering features of the Arctic: the backscatter from the 
trailing edge of an open Arctic lead and the backscatter 
from Arctic pressure ridges. 

The results of these experiments can be summarized as 
follows: 

1. The Arctic ice canopy and the laboratory model are 
acoustically non-rigid, since they do not have an 
infinite specific acoustic impedance mismatch with 
water. The consequence is that a Significant amount 
of acoustic energy incident upon such a surface will 
penetrate and be mode-converted to solid-borne 
vibrations within the solid material. 

2. Due to the effects of mode-conversion, diffraction 
theory alone is not sufficient to predict the amount 
of energy which will be backscattered from the 
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trailing edge of an open lead. Some of the acoustic 
energy incident on the underside of the Arctic ice 
canopy will generate solid-borne compressional and 
shear waves which propagate along the plate and 
reradiate energy back into the water upon reaching a 
discomt 1 nultey7 When superimposed with the 
diffraction, this reradiated energy results in 
backscattered signals larger than predicted by 
diffraction theory alone. Individual components of 
the backscattered signal were identified by 
establishing a source/receiver geometry which 
separated them by time of arrival. Additionally, the 
presence of head wave radiation from a smooth floating 
plate was discovered and measured. 

The backscatter from Arctic pressure ridges is 
affected by the amount of mode-conversion. Diffraction 
occurs from the corners of the ridge, but a signifi- 
cant amount of the incident acoustic energy also 
produces mode-conversion induced reverberation within 
the ridge which reradiates energy. The majority of 
the backscatter occurs in the least time region of the 
ridge. Also, for ridges large compared to an acoustic 
wavelength, reflection and the reradiation of ridge 
reverberation dominate the diffracted portion of the 


backscatter. Consequently, the characteristic 
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behavior of the backscatter resembles the directivity 
pattern for rectangular piston radiation from the 


slant face of the ridge that is ensonified. 
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APPENDIX A: PLANAR TRANSDUCER FOR MODE CONVERSION STUDY 


When the scattering characteristics of sound from the 
acoustically significant features of the acrylic Arctic ice 
model are well understood, many of the same experiments will 
be repeated in a layered waveguide to study the propagation 
and scattering of sound in a shallow water "Arctic" 
environment. In preparation for those experiments, it was 
necessary to build a planar transducer capable of effective 
transmission of normal modes of vibration in a waveguide. 
There was also an application for this transducer in the 
mode conversion study experiment. 

A layered waveguide is a water layer which is bounded on 
the bottom by a half space of sediment, and by a half space 
of air on top, with the characteristics described by Clay 
and Medwin (1977). Effective transmission of a single mode 
in a waveguide is possible if the driving elements are 
located to reinforce the mode pressure maxima, and if the 
phase between adjacent maxima is shifted 180°. Typical 
vertical pressure profiles for the first three normal modes 
in a layered waveguide are shown in Figure A.l. Since a 
greater number of transducer elements will provide better 
control of the pressure profile, a twelve channel amplifier 
was designed to provide twelve independent channels with 


separate phase and amplitude control from a single input. 
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lege Alls Vertical pressure profiles, Z,(z), for the 
First three normal modes ina layered 
waveguide. (m is the mode number) 
The term "phase" is loosely applied, since only a 180° phase 
shift 1s required between elements to produce the necessary 
excitation. Figure A.2 is an electrical schematic of the 
amplifier design. The amplifier uses 324C general purpose 
operational amplifiers (with four opamps per chip). The 
merst Chip uses only two of its opamps, and provides the 
Signal to the inverted and non-inverted buses of the 
amplifier. The 100K potentiometer is used to ensure the 
voltage on each electrical bus is identical. The second 
stage amplifiers provide amplitude control via the 250K 
potentiometer, and phase control via the toggle switch for 
each channel. 
Construction of the transducer was made possible by 
materials provided by the Raytheon Corporation. The active 
material for each element is the semicrystalline polymer 
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polyvinylidene fluoride (PVDF). MEwasmechosen fOr 1ts 
flexibility and toughness, its good acoustic impedance match 
with water, and for the ease with which it can be cut and 
handled. Each element was constructed by gluing two strips 
(3/8" x 24" x 0.022") of the copper-plated PVDF material to 
a 3/8" brass bar to increase the volume displacement, and 
therefore the projecting power, of each element. The twelve 
elements were stacked and glued in an aluminum casing with 
cork strips positioned to physically isolate elements from 
the caSing and from each other to damp vibrations which 
might be transmitted between elements. Finally, the entire 
transducer was sealed with Emerson and Cummings Eccothane 
7534, a low viscosity, medium hardness encapsulant to 
waterproof the transducer. A schematic of the transducer 
design is shown in Figure A.3. 

The performance of the transducer was tested by making 
near-field pressure response meaSurements in air, applying 
the theory of complex apertures described by Ziomek (1985). 
The complex frequency response in the near-field was 
measured by positioning a 1/2" Bruel and Kjaer condenser 
microphone at a height of 1/2" over 276 uniformly spaced 
positions across the active face of the transducer. The 
measurements were repeated for three different frequencies 
of interest and with various amplifier configurations. The 


transducer response at 20 kHz for a mode 1 and mode 2 
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excitation is shown in Figures A.4 and A.5. By the theory 
of complex apertures, if the frequency response can be 
measured for a given electrical input, then the acoustic 
response of the transducer is known. 

The PVDF transducer was used for an application in the 
mode conversion study experiment. With all of the elements 
in phase and equally amplitude weighted, the transducer was 
positioned below the smooth acrylic plate at the critical 
angle for compressional wave excitation in the plate. With 
compressional waves efficiently excited by the plane waves 
of the source, it was possible to effectively measure the 
plate compressional wave reradiation into the water from the 
end of the plate and as head waves from the water-acrylic 


materface. 
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APPENDIX B: FLEXURAL WAVE THEORY 


Flexural waves in a plate propagate by a bending motion 
about the neutral, undistorted axis. Portions of the plate 
above and below the neutral axis undergo a combination of 
compressional and shear deformation as the wave propagates. 
Flexural waves are also known as bending waves. Different 
frequencies will travel at different speeds, so propagation 
of flexural waves in a plate is dispersive, and a waveform 
will change its shape with distance of travel. Flexural 
waves in a plate will be generated for frequencies kh < l, 
where k = horizontal wavenumber and h = plate thickness. 

According to Ross (1976), when classical elastic theory 
1s applied to solid plates, longitudinal waves have a phase 


speed given by: 
Y As 


a 0 (1 =mie7) 

The plate longitudinal phase speed depends only on 
properties of the material. Vibrations in a thick plate are 
guite complicated, but if frequencies are restricted to the 
small kh criterion, then the analysis of thin plate 
vibrations is possible. For a harmonic disturbance in the 


plate, the Timoshenko equation for beams is used to describe 


the flexural motion: 


SZ 


Z Zee 2 2 2 4 
dt 7 1 12 at? 12 Cp * at” 








where h is the plate thickness (meters), c is the 
longitudinal or compressional wave speed in the plate, y 1s 
the plate displacement perpendicular to the surface, and lp 


is the shear parameter for the plate. The shear parameter, 





is a non-dimensional number which relates Young's Modulus 
Mme to the Shear Modulus (G). Phemconstant so f 
proportionality for a plate, Koc. 76 (1 + 92750), accounts 
for physical distortions in the plate due to the effects of 
shear. There are four solutions to the Timoshenko equation, 
two of which are imaginary, and the general solution can be 


expressed as (Ross 1976): 


v= ent eee Coie ae Dewy ene 
where w = 21f, k is the flexural wave number w/v¢, v¢ is the 
flexural wave phase speed, j =ik, i = ie and By C, 


and D are the complex coefficients found from the boundary 
Conditions. 

The flexural wave group and phase speeds for harmonic 
excitations of a thin plate can be found by substituting the 


general harmonic solution into the Timoshenko equation for 


3 


thine bceese The resulting equation, when rearranged to 


form a quadratic equation in terms of radian frequency, 


W, 
aS 

Py h? h? h? 

—_—_—— wt = Cee — k*? + 1 we + Cp” —j)k* = 0 
Cp” re 2 12 


Solving the quadratic equation for the flexural wave phase 


speed results in 


i 
ve? — - B + J B? ~ snc | (B.1) 


2A 
2 
l ph 
A = +e ; 
B = (fp, + 1)h?/12 
- = 2 2 
e = Cp hee 


It has been shown theoretically by Ross (1976) and 
experimentally by Denny and Johnson (1986), that in the low 
frequency limit, the flexural wave phase speed is 
proportional to the square root of the frequency: 

WCph *5 


vf ——— 
J12 
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Also, at high frequencies the flexural wave speed becomes 


constant, so flexural waves are then non-dispersive: 


2 


c 
Ve? = = = 0.92 cg 
ie 
where Cg = the shear wave speed in the plate. Siento kae dy ; 


when Equation B.1l is solved for the flexural wave group 
speed, it 1s seen that the group speed is twice the phase 


speed at low frequencies: 


wCph 5 


Jia 
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and that group and phase speeds are equal at high 


frequencies, 
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